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Abstract: 

In this contribution we show sufficient conditions for simultaneous unique identification of unknown 
spacewise coefficients and heat source in a parabohc partial differential equation given additional final time 
measurements. Our approach is based on density, in suitable spaces, of the corresponding adjoint problem. 

A second issue of this paper is the regularization approach. The sequence of approximated solution 
is obtained by coupling the nonlinear Landweber iteration with iterated Tikhonov regularization. We 
show that the parameter-to-solution map satisfies sufficient conditions to prove stability and convergence 
of approximated solutions for the identification problem. We use a unified discrepancy principle as the 
stopping criteria. 

Finally, we apply the developed theory in the inverse identification problem of unknown parameters 
(perfusion coefficient, metabolic heat source) for the identification of tumor regions by thermography. 

Keywords: uniqueness, thermophysical parameters and source identification, iterative regularization, 
parabolic type equation, final time measurements. 

1 Introduction 

Coefficients identification inverse problems have the characteristic of being ill-posed [4]. In other 
words, typically, solutions for such problems may fail to exist, may not be unique or be unstable 
under errors in the input data. The issue of existence can be relaxed by considering generalized 
solutions. On the other hand, uniqueness and stability are crucial for obtaining a reasonable solution 
for the coefficient identification inverse problem, theoretically as well as in terms of numerical 
approximations. 

Uniqueness in inverse problems have been studied for a long time. In particular, for coefficient 
identification in parabolic equations, see [26, 8] and references therein. However, only recently some 
attention has been given to the uniqueness of coefficient identification in parabolic type equation 
with final time measurements [2, 5, 8, 26]. A method that became very popular recently to prove 
uniqueness and conditional stability in coefficient identification inverse problems are Carleman type 
estimates, e.g. [8, 13, 26] and references therein. In particular in [26] there is a complete overview 
of Carleman estimates for parameter identification of inverse parabolic problems. Recently, in [2], 
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the authors proposed using Carleman type estimates for identification of spacewise source and heat 
conductivity in a heat type equation for a given additional final time measurement. For the one- 
dimensional heat equation, in [2], density properties of the associated adjoint problem, in suitable 
spaces and an additional assumption on the final measurements was used to guarantee uniqueness 
identification of spacewise source and heat conductivity. The duality type method was used before 
in [6, 5] for proving uniqueness identification of spacewise coefficients in parabolic type equations. 
However, in [6, 5] the author does not consider simultaneous identification of coefficients and heat 
source. 

Likewise uniqueness, regularization approaches for parameter identification in parabolic partial 
differential equations also have a long history and a full overview becomes almost impossible. See, 
for example, [15, 9, 4, 8] and references therein. We remark that, in general, parameter identifi- 
cation inverse problems are nonlinear, even if the forward problem is linear. Therefore, to prove 
convergence and stability of iterative regularization methods in this context, ones need to prove 
some nonlinearities conditions for the parameter-to-solution map [12], that are, in general, hard to 
be verified in practice. 

Summarizing, the main contributions of this paper are: We first use the density of the 
associated adjoint problem to prove sufficient conditions for the simultaneous uniqueness identifi- 
cation of spacewise heat source and coefficients that multiply the lower order term (or lower-order 
derivative, see Remark 2) of a parabolic type equation, with some assumption in the given final 
data (see equation (3) below). It is worth noting that we do not have restrictions of space dimen- 
sion, except that necessaries to prove existence and uniqueness of a solution of partial differential 
equation. Such restrictions are related to the smoothness of the initial and boundary conditions 
and the smoothness of coefficients of the partial differential equation. 

We also propose a iterative regularization method that consists in to couple Landweber and 
iterated Tikhonov regularization strategies. We prove properties of the parameter-to-solution maps 
(see definitions (6) - (8)) that are sufficient conditions to show convergence and stability of regu- 
larized solutions with respect to the noise level. We use a stopping criteria given by a unified 
discrepancy principle. Such discrepancy principle has the characteristic of reduce computational 
effort (See Remark 3). 

Finally, we apply the theory developed before to prove uniqueness and provide and regularization 
approach for the identification of the blood perfusion rate and the metabolic heat generation in a 
thermography application for melanoma diagnoses. 

The paper is organized as follows: In the remaining part of this section we introduce some 
notations. In Section 2, we present the model problem and define the parameter-to-solution map 
associated with the identification problem. In Section 3, we prove the uniqueness identification of the 
spacewise pair of coefficient and source from the additional final time measurement in the parabolic 
partial differential equation model. In Section 4, we prove properties of the parameter-to-solution 
map that guarantees the convergence and stability of the iterative regularized solutions, that will 
be proposed in Section 5. In Section 6, we present an application of the theory developed early in 
a melanoma diagnoses from thermography. Finally, in Section 7, we present some conclusions and 
future works. 



2 



Notations: By L^{Q) for 1 < p < oo, we denote the usual space of p-integrable functions on Q 
witli tlie usual norm IHI^p(q)- The space L°°{Q) is the standard L°°-space. We denote by W^'^^fl) 
the standard Sobolev space on Q with generalized derivatives of order < k in L^(f2). In particular, 
for p = 2 we have the Hilbert spaces H^{Q). 

Let T > be fixed. Define the measurable function u{-,t) : (0, T) — y X, where X is a Banach 
space. We denote by C([0, T]; X) the space of continuous mappings u{-,t) with the usual norm and 
by L^{{0, T); X) the space of functions such that 

T 

\\u{-, < oo . 



2 Model Problem 

In this paper, we are considering a thermal-physical model in a non-homogeneous and non-isotropic 
body, denoted by fl occupying an open, bounded and smooth domain in M", described by the 
parabolic type partial differential equation 

Ut — L{a, b, c)u = f{x) in f2 X (0, T) 

t) = for (x, t) Ednx (0, T) (1) 
u{x, 0) = f{x) for X G , 

for a time interval (0, T) with T > 0, where 

L(a, b, c)u = V ■ (a(x) Vm) — b{x) ■ Vu — c{x)u , (2) 

is a linear elliptic differential operator of second order with all the coefficients time independent. 
Moreover, a, c are strictly positive real valued function in L°°(f2) with < a < a{x) and < 
c < c{x), for X G 6 is a real valued vector function sufficiently smooth. / G L^(f2) is the 
spacewise heat source. For simplicity, we assume that the given initial temperature distribution 
(/? G H'^{Q) n Hq{Q), and that 6 = 0. Hence, we can assume that information that u is identically 
zero on the boundary of Q. We will denote L{a, b, c) = L{a, c). 

Moreover, we assume given the additional final temperature measurement g G H'^{Q) fl H^{Q) 
satisfying 

u{x,T) = g{x) foT X eft, T>0. (3) 

Since the parameters in (1) are sufficient smooth, we can, formally, define the adjoint of the 
partial differential operator L{a,c) as 

L(a, cYv = V ■ {a{x)Vv) + V{b{x)v) - c{x)v . (4) 

Given the assumptions on the parameters in the model (1), follows from the Hille-Yosida 
Theorem [27] that the operator —L generate a strictly dissipative and contraction Cq semigroup 
{G{t)}teR+ in L^in) with V{L) = {u : u e H^iil) , Lu G L\n)}. Hence, ||G(T)|| < 1. Note that, 
since Hq{Q) is compact embedded in G{t) is a compact operator for every t > 0. Moreover, 
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it follows from classical results on parabolic partial differential equations, e.g. [16, 14] that there 
exists a unique solution u G C^{{0,T), H^{n) n H^i^)) of (l)-(3) with 




dt < oo. 



(5) 



Remark 1. i) Is possible to extend the result to a weaker assumption on the data, but one then 
has to consider appropriate weak formulation of (l)-(3), [14]- 

a) Since we have time-independent coefficients and source, it follows that the solution u to (1) 
with additional data (3) is analytic in time. This means that u has derivatives of all orders 
with respect to t [8, 25] . Moreover u is at least continuous with respect to time, and, therefore, 
we can conclude that pointwise evaluation in time makes sense. 

The inverse problem that we are interested here is recover the pair of spacewise parameter and 
source {c{x),f{x)) in (1) from the additional final time measurement (3). 

Assuming that b = b{x) is not zero in (1), one can recover the pair of spacewise parameter and 
source {b{x), f{x)) in (1) from the additional final time measurement (3). However, the uniqueness 
result, in this case, follows from a very similar argument (See Remark 2) developed with this 
contribution and we will not fix in this problem here. 

2.1 The parameter-to-solution map 

In this subsection, we introduce the parameter-to-solution map associated with the parameter iden- 
tification problem discussed previously. For now, we will consider the following admissible set of 
spacewise coefficient and heat source: 

Definition 1. We denote the admissible set of functions as 



Moreover we denote by T>c{¥) C L°°{Q) and 'DfiF) C the restriction ofV(¥) to the first and 

second component of the pair (c, /), respectively. 

Note that, since Q is bounded, V{¥) is a convex and closed subset of x L^[Q). However, 

since L°°{Q) can not be continuously embedding in H^{Q) for ^2 C M^, V{¥) has no interior points 
when equipped with the H^{Q) norm. This will not affect the convergence analysis that follows. 
See Section 4 for details. 

Moreover, let u e C([0,T]; H^{n) D H'^{n)) be the unique solution of (1), with (c, /) E V{¥). 
Then, it follows that the restriction u{x, to) is well-defined for < to ^ ^- Therefore, the restriction 
of u to r2 X {T} exists, it is, u{x,T) = g{x) is well defined and, moreover, the nonlinear operators 



I?(F) := {(c, /) G X L\n) s.t. < c < c(x) <c,a.e. in Q} . 




¥{c,f) = g{x) 



(6) 



is well-defined. 
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A second operator equation that we are interested in is the restriction of the operator F defined 
in (6) in Vf{¥). Indeed, it introduces a family of operators, parameterized by c G V(¥) defined by 



/^F,(/):=F|^^^,,(c,/) = ^?(x). (7) 

Note that, for any fixed c G P/(F), the operator Fc is nonhnear, unless ip = 0. 

Finally, we introduce the restriction of the operator F to define a family of nonlinear 

operators, parameterized by / G defined by 

Af -.v^iw) c L'^in) — > L'^{n) 

c^Aj{c):=¥\^^^^^{cJ) = g{x). (8) 

In practical applications the final temperature (3) is, in general, not known exactly. One is given 
only approximate measured data G L'^{^), corrupted by a noise level 6 > satisfying 

Therefore, given the measurements g^, the inverse problem is: Find (c(x), /(x)) G V(¥) such that 

F(c,/)=/, for / satisfying (9) . (10) 

The inverse problem is ill-posed in the sense of Hadamard [4]. In other words the solution of the 
inverse problem is unstable with respect to noise data, it is, small perturbation in the given data 
implies in large perturbation on the parameter space. Examples of instability of reconstruction the 
coefficient c can be obtained in [6]. For instabilities examples in the reconstruction of / see [10, 22]. 
Hence, beyond uniqueness, some regularization methods need to be used to guarantee stability of 
the parameter and source reconstructions, given a set of noisy data. 



3 Uniqueness for the spacewise coefficient and source 

In this section, we wish to prove uniqueness of the spacewise coefficient c and source / in (1) by 
additional measurement (3). We will use a similar approach of [2, Section 3] for the identification 
of the heat conductivity and heat source for the one-dimensional version of the heat equation (1). 
The derivation of the uniqueness result is based on a completely different technique than Carleman 
estimates [2, 8, 26]. Indeed, the technique is based on results of density, in certain function spaces, of 
solutions of the corresponding adjoint problem [6, 5] and the unique continuation principle [16, 14]. 
Moreover, the proposed approach is different of the maximum principle used in [8, Section 9.1]. It is 
worth remark that, differently of [8, 6, 26], we wish to show uniqueness of the spacewise coefficient 
c and source /, simultaneous, from (l)-(3). 

The steps for proving uniqueness of the identification of the pair of parameters {c{x), f{x)} for 
giving initial and final data in (l)-(3) are outlined as follows: 
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Preliminary results: Denote by m = m(ci,/i) and v = u{c2, f2) the respective solutions of (1) 
with additional data (3), for coefficient and source satisfying the Definition 1. Then, for linearity 
of (1) the difference w = u — v satisfies 

wt - L{a, c,)w = (ci(x) - C2{x))v + - /2(x)) inQx (0, T) (11) 

with homogeneous initial, boundary and final conditions. 
Let we invoke the adjoint problem of (11), that reads as 

ipt + L*{a, ci)ip = in 1] X (0, T) 

i;{x, t) = for {x, t) e dn x (0, T) (12) 
ip{x,0) = for X e 1], 
ifj^x^T) = fi{x) for X G f2, , 

where fi{x) is an arbitrary function in Cq{Q). 

For properties of solutions of the adjoint equation (12) we have: 

Lemma 1. Let the Assumption on the coefficients, source, initial, boundary and final conditions in 
(l)-(3) hold. Then: 

i) For any function fi{x) G Cq{Q), there exists a unique solution ip{x,t;v) G C^{{0,T); Hq{Q) fl 
H'm of (12). 

a) For any function fi{x) G Cq{Q), the following relation holds 

ip{x,t; ^{x))F{x,t)dxdt = 0, (13) 



'0 Jn 

where w is a solution to (11) with right-hand side F{x,t) := (ci(x) — C2(x))f + (/i(x) — /2(x)). 

Hi) For fi{x) ranging over the space Cq{Q), the corresponding range of ilj{x,t; fi{x))\t=r is every- 
where dense in L^(f2) at any time t = r, < t < T . 



iv) Given that 



T 



■ip{x, t] f^{x))^{x, t)dxdt = 



'0 Jn 

for fi{x) ranging over the space Cq{Q), then $(x, T) = , a.e. x G f2 . 

Proof. Item i) is a well-known result for parabolic equations, see, for example, [16, 14]. 

Item ii) follows immediately by multiplication of (15) by ip and integration by parts. 

Item iii) and Item iv) are consequences of [5, Lemma 2-3, pg 318] applied for a multi-dimensional 

case. □ 

Lemma 2. Let w satisfying (11). Then w{x,t) = almost everywhere for {x,t) G x [0,T]. 
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Proof. Note that, since w has homogeneous boundary conditions, it follows for standard parabolic 
theory that w{-,t) G Hq{Q), for any t G (0,T). From Remark 1 ii) we have that w has derivatives 
of all orders with respect to t and time pointwise evaluation makes sense. 
Applying Lemma 1 Item ii) in combination with iv) in (11) we get 

[(ci(x)-C2(x))t;(x,T) + (/i(x) + /2(a:))] = 0, a.e. xeQ. (14) 

Using this in the first equation in (11) in combination with w{x, T) = for a.e. x in Q, we conclude 
that Wt{x,T) = 0, a.e. x in Q. 

Define z = Wt- Since coefficients and source are time independent, we have that z satisfies 

Zt — L{a, b, cijz = (ci(x) — C2{x))vt in Q x (0, T) 

z{x, t) =0 for (x, t) eOnx (0, T) (15) 
z{x, 0) = 9i{x) for X G , 
z{x, T) = for X G Q. 

Splitting this problem into two, one with zero right-hand side and with initial condition 6i and one 
with the given right-hand side and zero initial condition, following the proof of [2, Lemma 2] one 
can conclude that the solution to the first one is identically zero, i.e. 6i{x) = 0. Therefore, since 
z satisfies a problem of the same kind as w we can again apply Lemma 1 Item ii) in combination 
with iv) in (15) to conclude that 

{ci{x) — C2{x))vt{x,T) = , a.e. x E Q . (16) 

Using this in the first equation in (15) in combination with z{x, T) = for a.e. x in Q, we conclude 
that zt{x,T) = 0, i.e. wtt{x,T) = 0, a.e. x in Q. Continuing this, putting zi = zt and deriving the 
problem for zi and applying the similar reasoning, i.e. Lemma 1 Item ii) in combination with iv), 
we find that Wtu{x, T) = 0. Further continuing this it is possible to prove that (9f^^^;)(x,T) = for 
/c = 0, 1, 2, . . .. From this and strong unique continuation results for parabolic equations [16, 14], 
we conclude that w{x,t) = for a.e. {x,t) G f2 x [0,T]. □ 



The uniqueness proof: We now have the required results in order to prove the main step in the 
uniqueness of (c(x),/(x)) in (1), with additional final data (3). 

Theorem 3. Let the Assumption on this paper holds. Moreover, assume that g{x) ^ fix) for 
a.e X G fi. Then the inverse problem (l)-(3) has a unique solution {m,c, /} with the coefficient 
c G L°°{Q), the heat source f G L'^{Q), and temperature u, with < oo. 

Proof. Follows from Lemma 2 that w is identically zero in f2 x [0,T]. This in particular implies 
that z = QuiVtx [0,T]. 

From the first equation in (15) we then have 

(ci(x) - C2(x))ft(x, t) = , for a.e. (x, t) eVLx [0, T] . 

Since the coefficients are independent of time, we integrating with respect to time, form to T and 
use the fundamental theorem of calculus, to get 

(ci(x) — C2(x))(5'(x) — f{x)) = (ci(x) — C2(x))(f (x, T) — v{x, 0)) = , for a.e. x G f2 . 
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From assumptions on g and if and we can conclude that Ci(x) = C2(x) also for a.e x G fi. 
Moreover, since w = (from Lemma 2) and ci(x) = C2(x), we have from (11) that 

/i(x) — f2{,x) = , for a.e. x eVL . 

The uniqueness of m, with follows from the standard theory of solution of parabolic partial 
differential equations [14]. □ 

It is worth to note that the argument in the proof of Theorem 3 goes beyond the proof of 
uniqueness in [2]. The main reason is that, for the lower order terms we do not have the influence 
of the divergent operator. Therefore the proof is still true in multidimensional heat equations. 

Remark 2. With similar argumentation of Theorem (3) is possible to prove uniqueness of {u, b, /} 
for the model (l)-(3). Indeed, some small modifications in the derivations of the steps above are 
necessary. The main difference is that in Theorem 3 we need the assumption that \ V g{x) — V '^{x)\ > 
for a.e X E fl. 



4 Properties of the Parameter-to-Solution Map 

Before introducing the iterative regularization in Section 5, we need to prove some properties of 
parameter-to-solution map defined before allowing us to obtain convergence, stability and regular- 
ization properties of approximated solutions. 



4.1 Continuity 

The first result in this direction is the continuity of operators F, A/ and Fc defined in (6) - (7), 
respectively. 

Theorem 4. The operator ¥ : 'D(F) C L^(f2) x L^(f2) — > -L^(i7), defined in (6) is continuous. 

Proof. Let {(c„, /„)} be a sequence in V(¥) converging to some (cq, /o) € V(¥) w.r.t. L^(f2) xL^{Q)- 
norm. Denote by u„ = u(c„, /„) and v = u{co, /o), respectively, the solutions of (1). As before, the 
difference w := Un — v satisfies 

Wt-V ■ (a(x)Vw) + Cn{x)w = (Co - Cn)v + (/n - /o) 

with homogeneous, initial, boundary and final conditions. 

Since v G L^{0, T; H^in)) C([0, T]; H^in)), for each t G [0, T] we have 

WtW — V ■ {anVw)w + Cn{x)wwdx = / {Cq — Cn)vwdx + / {fn — fofWdx 

fi Jn Jn 

The Green's formula [16], implies that 

^^lkWllL(o) +a||Vu'(t)||j2(t^) + c||u7(t)||5^2(^) - \cn- co\\v{t)\\w{t)\dx + Jjfn- fo\\w{t)\dx, 
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where we used the homogenous initial and boundary conditions of the equation satisfied by w. 
Moreover, from Theorem A. 2, there exists some e > such that Ibl^)!! < C, for q = 2 + e. 

Since the apphcation t ||f (^) || continuous {v G C{[0,T],W^'''{Q))), we have that is 
uniformly bounded for t G [0,T]. 

Let + = 2~^. Using the Holder inequality with + q^^ + 2^^ = 1, we have 

<cf||c„-co||^,(f^) + ||/„-/o||i2(f^)) \\w{t)\\^,^^y 



Since, for each t G (0,T) G Hq{Q), it follows from the Poincare inequality that ||w;||^2(q) < 
Cill Vuijl ^2('j^')- This, together with the Young inequality with i [16], yields 

d /.X ,,2 ,11 /,M|2 



rmn{2 \Cia,c} { ^lk(t)IU2(t^) + \\w{t)\\ L^i^n) 



< J (||C„ - Co||^,(^) + Wfn - /o|L2(f,)) + C5||u'(t)||'2(^) 



Let e > such that Ce > min{2 \Cia, c}. Using the Gronwall inequality, it follows that 

Ik(t)|li2(f^) <e(||cn-Co||^p(f^) + ||/„-/o||^2(f^)) e'^* 

for all t G [0,T]. Given the continuity of the solution of (1) with respect to t the inequality holds 
for t = T. 

Now the conclusion follows form Lemma A.l in Appendix. □ 

As a corollary of Theorem 4 we conclude the continuity of operators Af and Fc, for any fixed / 
and c in T>(¥), respectively. 

Corollary 5. For each fixed c G T>{¥), the operator Fc defined in (7) is continuous in L'^{Q). For 
each fixed f G ViF), the operator Af defined in (8) is continuous in 



4.2 Frechet Derivative and Tangential Cone Condition 



An important result to guarantee convergence of iterative regularization methods for nonlinear 
inverse problems is the local tangential cone condition. We will prove such properties in the next 
two propositions. 

Proposition 6. For each fixed c G "DclW) the operator F^ is Frechet differentiahle. The Frechet 
derivative is Lipschitz continuous and satisfies the local tangential cone condition. In other words, 
for each f,fE Bp{fo) C Df(¥), there exists a < ri < 1, such that 



F^Cf) - F,U) - KU)U - f) 



L2(Q) 



FcU) - F,{f) 



L2(Q) 



(17) 
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Proof. Let h G By linearity of (1), the sensitivity u' ■ h = u{a, f + h) — u{f) satisfies 

iu' ■h)t-W ■ {aViu' ■ h)) + c{x)u ■h = h, (18) 

with homogeneous initial, final and boundary condition. Follows from the standard parabolic theory 
that there exists a unique solution C{[0,T], Hq{Q) x H'^{Q)) of (18) and 

11^ ■ ^\\c([o,T],m{n)) - • 

Let {G{t)}t^^+ be the semigroup generated by the differential operator —L{a,c). Define the 
linear operator K{t) : L\n) — ^ L\{0,T); H^{n)) n C([0,T]; H^{n)) by K{t)f = f^G{t - s)fds. 
Note that the solution u of (1) is formally given by u{x,t) = G{t)g{x) + K{t)f{x). Therefore, 
solving the operator equation (7) is equivalent to solving 

K{T)f = u{c,x,T)~G{T)g{x). 

In other words, {,Fl.{f)){h) = K(T)h = u'{x,T) ■ h, where u' ■ h is the unique solution of (18). 
Therefore, follows from the linearity and continuity of K that is Lipschitz continuous and satisfies 
the tangential cone condition. □ 

Proposition 7. For each fixed f G Vf[¥), the operator Af is differentiahle in the direction k, such 
that c + K G Dc(¥). The derivative can be continuously extended as a linear operator to H^{Q). 
The extension is Lipschitz continuous. Moreover, the local tangential cone condition is satisfied. In 
other words, there exists p > and < rj < 1 such that for each c,c & -Bp(co) C Dc{¥), 

\\Aj{~c) - Af{c) - 4(c)(g - c)||^,^^^ < T^WAfid) - Af{c)\\^,^^^ . (19) 

Proof. Let k G Dc{¥). Then, by linearity and continuity with respect to the coefficients of solutions 
of equation (1) we have that the directional derivative ■ (k) in the direction k such that c + k G 
V{¥) satisfies 

■ ii^))t - L{a, c)m'(,) ■ (k) = KM , (20) 

with homogeneous initial and boundary conditions. It follows from standard parabolic partial 
differential equation theory [16, 14] that there exists a unique solution uj^^ ■ (k) G C^((0, T); Hq{Q) fl 
H'^{Q))nG^{[0, T]; L^{Q)) of (20). Hence, the u[^^ ■ {K){t) make sense, for every t G [0, T]. Moreover, 

hU ■ (^)|lcO([0,T],L2(n)) - ll'^"llL2(nx(0,T))- 

Since k is time-independent, follows from the Cauchy-Schwarz inequality that 

hU ■ ('^)(^)IL2(f,) <CiT)Mm{n)\Mi- 

Therefore, the directional derivative u|^^-(k) can be extended to H^{Q) as a bounded linear operator. 

To prove the Lipschitz continuity, let c, c G P(F) and u = u{c), u = u{c) the respective solution 
of equation (1). Then the difference w = ■ (k) — u'^^-^ ■ (k) satisfies 

Wt — L{a, c)w = (c — c)u + k{u — u) . 
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As before, we have 

11^2(0) < l|c — ^ + — ^t|li2((o,T)xQ) ■ 

With the same argumentations as in Theorem 4, we have \\u — '^||j;^2((oT)xn) — ^W'^ ~ ^Wn'^in) ^^'^ 
the Lipschitz continuity follows. 

Moreover, from the linearity of equations (1) and (20) we get that v = u{c) —u{c) — u'{c) ■ (c — c) 
satisfies 

Vt — L{a, c)v = (c — c){u{c) — u{c)) , (21) 

with homogeneous boundary and initial conditions. Using similar argumentation as in the proof of 
Theorem 4, we obtain that 

\ d f 

2^IK'Wlli2(n) + «l|Vt'(^)lli2(f,) + c|I^WIli2(f,) < / {c-d){u{c)-u{5))v{t)dx 

11^ ^IlL'^fn) ~ 2 2 

< 1 II^(C) -^^(5)||^2(f^) +e||f(t)IU2(f^), 

where we use the Young inequality with e. 
Let we take e < c/2. Therefore, we have 

WvmL^^n) < C\Hc) - w(5)IU2(f,) , Vt e (0,T) 

where C = C{\\c — c)|| and some constants that are independent of the solution of (1). Denoting 
?7 := C(||c — c|| ), we have that 



Given the continuity of v and u with respect to t, the inequality holds for t = T. □ 

Adjoint of the Prechet derivative: Let we finish this section making the calculation of the 
adjoint of the Frechet derivative of the operator defined in (7). 

Lemma 8. Let Vf G L'^{Q). Then the adjoint of the Frechet derivative F^{f) denoted by {F^{f))* : 
L^(f2) — L^(r2) is given by 

iKif)nrf) = -V{x,0), (22) 

where V G C{[0,T], Hq{Q)) is the unique solution of 

Vt + V- {aVV) + cV = rf, m (0, T) x , (23) 
V{x,t) = 0, (x,t) G 9fi X (0,T) 
V{x,T) = 0, for X eft. 

Proof. Existence, uniqueness and regularity of the solution of (23) follows from standard parabolic 
theory [14]. Note that equation (23) is the adjoint equation of (18). Therefore, the assertion follows 
directly from integration by parts. See the details of the calculations in [10]. □ 
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5 An Iterative Regularization Method 



In this section, we propose an iterative regularization method to regularize the solution of the 
inverse problem (6). It consist in a coupled Landweber - iterated Tikhonov regularization approach 
given by the iteration 

Given Co = 4, fo = fo e V{¥) , for A: = 0, ■ ■ ■ , A:, 

fLi = ft + liF'AftWig' - F^iUt)) (24) 

A: ^ 

2 



4j^i G argmin Ja{c 



II <5II^ 



where, is the stopping index, determine by the stopping criterion using the discrepancy principle 
and the relaxation parameter r is such that 



|F(cf+i,/,V) -/||^,^^^ <t6< ||F(ci„/,^) -/||^,^^^, (25) 



r > 2i±^ . (26) 
1 — r] 

Note that, if we have noise free data then /c* = +oo. In this case we drop the index 6 in (24). 
Moreover, in (24) the positive parameter 7 is a scaling factor to enforce convergence of the 
Landweber iteration [12, 4]. As a consequence of Lemma 8, we have 

fLi = ft + lV,{x, T) = fi + 7i^*(T)(/ - Ukix, T)) , 

where Vk{x,t) is the unique solution of (23) with rj = Uk{x,T) — and Uk{x,T) = u{cl, fl) is 
the unique solution of (1) with c = cl and / = /f, respectively. Therefore, is enough that 7 < 
||i^(T)||~^. Since the operator —L generate a contraction semigroup, it follows that ||/('(T)|| < T. 
Therefore, is enough that < 7 < T~^. This estimate in not sharp. 

It is worth noticing that in iterated Tikhonov approach, the parameter a do not play the rule 
of the regularization parameter [1]. In this case, we can choose any a > 6((5/p)^, where p is the 
radius (fixed) of the ball around cq. 

Remark 3. In this remark we will discuss some point about the proposed iteration. 

• The proposed algorithm (24) is a type of Kaczmarz strategy [11]. However, it is not the same 
Kaczmarz iteration proposed before in [7, 3, 1], since the unknown is a pair of parameter (c, /) 
and not a single parameter. 

• Moreover, the proposed iteration is not the same as using a Landweber iteration [12] for the 
coordinate fk and the iterated Tikhonov [1] for the coordinate Ck in the pair (cfe, fk), since the 
iteration Ck+i depends of the iteration fk+i o-s we can see in the iteration (24). 

• We can mix some other type of iterative regularization methods in order to regularize the pair 
of parameters {c,f). However, the choice iterated Tikhonov in the second line in (24) imply 
that we can use the uniform discrepancy principle. The advantage of this choice is that, in 
each iteration, we only need to evaluate one time the residual, indeed, at the end of the cycle in 
the algorithm (24) . See also (28) below. It saves significatively computational effort, compared 
with a discrepancy principle defined for each one of the lines of the system (24). 
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5.1 Convergence Analysis 

We start the analysis of the proposed algorithm with the following result that imply in the well 
posed of the iterative Tikhonov regularization. 

Lemma 9. For each f G Df(¥) fixed, there exists a minimizer of the Tiknonov functional Ja defined 
m (24). 

Proof. Note that -D/(F) is convex and closed in Li^iVt). Therefore it is weak sequentially closed. 
Now the proof follows immediately from the continuity oi Af given by Corollary 5. □ 

Given the iteration formula in (24), we conclude that 



(27) 



As usual for nonlinear Tikhonov type regularization, the global minimum for the Tikhonov func- 
tionals in (24) need not be unique. However, in [1] was proved that, for exact data, is possible to 
obtain convergence statements for any possible sequence of iterates, and we will accept any global 
solution. For noisy data, a (strong) semi- convergence result is obtained under the assumption that 
Af has a Lipischitz Frechet derivative as we have proved in Proposition 7. 
Given the minimality of c^_,_^ in the iteration (24), we have 



\¥(cl 



9 



L2(n) 



(28) 



Therefore, one important consequence of (28) is that, if the unified discrepancy principle (25) is 
not attained at the iteration k + 1, then the standard discrepancy principle for Landweber iteration 
also is not attained, it is, while 



r6 < \\¥{cU„ f,^,j 



9 



lL2(n) 



then t6 < 



^cdfk 



5 ^ 
k+1 J 



9 



L2(f7) 



(29) 



Because of this inequality, we call the discrepancy principle (25) a unified discrepancy principle. 

Now, we are able to prove the convergence and stability of the iterative regularization method 
in (24). 



Theorem 10. Let (cq, /o) 



-0' 



/o) G V(¥) and the operators, ¥, Fc and Af as defined in (6) - 



(8) and t as in (26). Then, for any {c*,f*) G T){¥) a solution of (6), the iteration given by (24) 
has the following properties: 

1. While \ W - ^^-^ 



L2(n) 



> t6, we have that 



||r 


- 1 

Jk+l\ 




< 


|r- 


fl 


lL2(f7) 


||c* 


- 1 


L2(n) 


< 


\c* - 




L2(n) 



(30) 



Moreover, if (cq, /o) G Bp{c\ /*) C V{¥), then {cl, fl) G B2p{c*, /*) for all k and 



fc,-i 



k^{T5r<j2y-F4f, 



k=0 
k=0 



< 



r\\f*-f^ 



0llL2(f1) 



L2(n) 



;i - 2r])r - 2(1 +r/) 



fc+i 



r c 



< 



'^0|lL2(n) 



[1 - 2r])r - 2(1 + 7/) ' 



VO < A; < A;* 



VO < A; < A;* 



(31) 
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In particular, if = g (i.e., 6 = 0), then 



I]ll5-^c,(/fc)|li2(t^) < oo and ^||c/-%^^(cfc+i)||^ <oo. (32) 



fc=0 k=0 



2. If there exist (c*, /*) G 5p((co, /o)), a solution of (6) and 6 = 0, then there exist a subsequence 
{ck, fk) given by (24) that converges to {c*,f*). 

3. In the noisy data case, if the iterations are stopped according to the discrepancy principle (25) 
and T is given by (26), then there exists a subsequence (that we denote by the same index) 
(^k{5,g^y fk{5,g^)) ^^^^ convcrgcs to a solution {c* , f*) of (6), as 5 0. 



Proof. Let cf. be fixed. Then, since F^s is continuous, compact (see Corollary 5) and satisfies the 
local tangential cone condition (see Proposition 6), follows from [12, Chapter 2] that the sequence 
fk ~ fki^k) given by the Landweber iteration satisfies the claim of the Theorem. 

Now, let fl be fixed. Then, since Ajs ^ is continuous, compact (see Corollary 5) and satisfies the 

local tangential cone condition (see Proposition 7), follows from [1] that the sequence = cl{fl) 
given by the iterative Tikhonov method satisfies the claim of the Theorem. 

Therefore the convergence, stability and regularization properties of the approximated sequence 
{of., fl) follows from a diagonal argument. □ 



6 Application in Thermography 

Nowadays is well known that the body surface temperature is controlled by the blood perfusion, 
local metabolism and the heat exchange between the skin and the environment. Changes in any 
of these parameters can induce variations of temperature and heat fiuxes at the skin surface. In 
particular, the apparition of malignant tumor imply in a highly vascularized skin region that lead 
an increase of local blood fiow. Consequently, in a local increases of the blood perfusion and of the 
capacity of metabolic heat source [24, 23, 17]. 

Modern diagnostics of melanoma location in the skin are non invasive. They use the skin 
surface temperature measurements. However, this technique requires the solution of inverse bio- 
heat transfer problem. This problem consists in the simultaneous identification of thermal and 
geometrical parameters of tumor. In applications, some of the parameters that are interesting are 
the perfusion coefficient and the capacity of metabolic heat source [24, 19, 17, 23]. 

From the mathematical point of view the heat transfer processes in the domain of biological 
tissue are described by the Pennes [24, 20, 28, 21] equations 

pCpUt - V ■ {aVU) - iOb{x)pbCb{Qo -U)= Qm{x) , in x (0, T) (33) 

where p, Cp, a denotes density, specific heat, and thermal conductivity of tissue; pb, Cb are density and 
specific heat of blood; Ub blood perfusion rate; Qm metabolic heat generation; Qo is the supplying 
arterial blood temperature and U the tissue temperature. Q is the body region around the melanoma 
location. Therefore, C M" for n = 2 or n = 3. For simplicity, we assume that the melanoma is 
located just below the skin and that we can consider Q C M^. For the case of f2 C the model is 
more complicated, principle, in terms of the boundary conditions [24, 20, 28, 21]. 
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Besides the thermal parameters and metabohc rate of tissue, the skin temperature is also deter- 
mined by many other factors such as the skin humidity, radiation emissivity of skin and parameters 
of surrounding air. These factors can be incorporated into the boundary condition at the skin 
surface. However, for simplicity, we assume that far from the tumor location, the heat effect 
of the tumor activity is insignificant. Therefore, the boundary conditions can be assumed to be 
constant and equal to Qq. Moreover, without loss of generality, let us consider the parameters 
p = Cp = pb = Cf, = 1. Denoting u = U — Qq we have that u satisfies 



where (p{x) imposes an initial spatial heating, c{x) = ujh{x) and f{x) = Qm{x), that we assume be 
smooth as in the above sections. We assume that the temperature measurement, at final time, is 
given as in the equation (3). 

Therefore, the theory developed before in this paper is applicable to the melanoma location in 
the body given the measurement on the skin surface in the following sense: 

Theorem 11. There exists a unique blood perfusion rate Ub{x) and a unique metabolic heat gener- 
ation Qm{x) for a given skin temperature measurement satisfying (3). 

Proof. Is a direct application of Theorem 3. □ 

Theorem 12. Given measurement of the temperature on the skin surface satisfying (9). There 
exists an iterative regularization method that generates a sequence of approximate solutions for the 
identification of the blood perfusion rate Uf,{x) and the metabolic heat generation Qm{x) that is 
convergent and stable w.r.t. noise in the data. 

Proof. Consider the iterative regularization approach given by (24). Then, Theorem 10 imply the 
assertion. □ 

7 Conclusions and Further Works 

In this work, we prove uniqueness of the spacewise parameter c and the source / in (1), for a given 
extra final time measurement (3). Moreover, we derive sufficient properties of the parameter-to- 
solution map to guarantee convergence and stability of approximated solutions obtained by the 
proposed iterative regularization method, if the stopping index is determined by the discrepancy 
principle (25). We also have analyzed the application of the theory developed here for a simplified 
version of the thermography model in melanoma identification. We are able to say that exists a 
unique blood perfusion rate and a unique metabolic heat generation in the simplified model (33) 
for a given final time measurement (3). 

The next step in this line is the numerical implementation. In particular the analysis of the 
three-dimensional Pennes equation was not totally covered. In the three-dimensional case there are 
also many numerical difficulties, beyond the theoretical, that needs attention [24]. 



Ut — V ■ (aVu) + c{x)u 
u{x, t) 
u{x, 0) 



fix), {x,t)enx{o,T) 
0, (x,t) G 9fi X (0,T) 
<^{x) , 



(34) 
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A Appendix 

In this appendix we will provide an auxiliary lemma that we need in the prove of Theorem 4. 

Lemma A.l. Let {ipk}keN be a sequence of functions, with c < ^Pk{x) < c for all x eVL. If ^ ip 
in L^{Q) for some p G [1, oo) , then ip^ ip in Lp{Q) for all p G [1, oo). 

Proof. We remark that, since Q is bounded, ipk G U\Q) for all 1 < p < oo. 

Assume — (/? in L'^{Q). 
Case 1.1 {2 < p < oo) 
For all n, / G N with n,l > 



Jn Jn Jn 



(pi\^dx . 



Hence, {ipk} is a Cauchy sequence in L^iVL). Therefore, — in ,2 < p < oo. 

Since, U'[Q) is continuously embedding in L^(f2) for 2 < p < oo, we have (p G L'^iVL) fl L^iVl) 
and 

llv^fe - ^l|L2{n) < C\\ipk - v\\LP{n) ■ 

By the uniqueness of the limit ip = (p. 
Case 1.2 (1 < p < 2) 

For all ?i, / G N with n,l > ko 

(fn- ^il^dx = / (\(pn- ^il'^j^dx < {meas{il)y'' i / \ipn-^i\'^dx 
In other words, {(fk} is a Cauchy sequence in Lp{Q) and hence, ip'' ^ ip in Lp{Q) for 1 < p < 2. 

Claim. <p G L'^{n)nLP{n). 

J \^\^dx<c(^j \ipk - 'fl'^ + l^kl^dx-^ < C (^J (^\Lpn - ipil^Y ^ j '''^'=1^^^^ 

< C (^eas{yL)f~l j \ipk - 'p\^dx + C-^ < oo . 

Hence, ip G L'^{VL). The continuous embedding of L'^iVt) in L^iVt), for 1 < p < 2, conclude the 
claim. 

Now, given e > 0, there exist /c G N large enough such that 

llv^ - <^||LP(n) < WVk - + ll'y^fc - < C\\(pk - <p\\mn) + llv'fc - vWLPin) < £ ■ 

Therefore, ip = ip. 

The arguments used in the proof of the reciprocal are similar to those presented above. Thus 
we will omit the proof. □ 
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The next theorem is a version of Meyers's Theorem [18, Theorem 1] adapted to our case. 

Theorem A. 2. [Meyers] Let Q and the coefficients {c, f) G V(¥) as in Definition 1. Then, there 
exists a po > 2 such that the unique solution u = u{a, f) of (1) belongs to L'^{0,T]W^'P{fl)) fl 
L\0,T;H',{n)) for anype[2,po[. 

Proof. It follows from the classical parabohc partial differential theory that that u G L^{0, T; ifg (r2))n 
C([0, T]; i^Q (fi)), [14]. Therefore, for each t g]0,T[, u{-,t) satisfies the following elliptic equation 

— V(a(x) Vn(a;, t)) + c(x)n(x, t)t = f{x) , x G ^2 

t) = on dVt . 

Follows from Meyers's Theorem [18, Theorem 1], that there exists a po > 2 such that u{-,t) G 
iy^'P(r2) for all p G [2,po[- It proves the assertion. □ 
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